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Many sedimentary rocks are formed by migration of sand ripples. Thin layers of coarse and fine 
sand are present in these rocks, and understanding how layers in sandstone are created has been a 
longstanding question. Here, we propose a mechanism for the origin of the most common layered 
sedimentary structures such as inverse graded climbing ripple lamination and cross-stratification 
patterns. The mechanism involves a competition between three segregation processes: (i) size- 
segregation and (ii) shape-segregation during transport and rolling, and (Hi) size segregation due 
to different hopping lengths of the small and large grains. We develop a discrete model of grain 
dynamics which incorporates the coupling between moving grains and the static sand surface, as 
well as the different properties of grains, such as size and roughness, in order to test the plausibility 
of this physical mechanism. 

PACS: 81.05.Rm, 91.65.Ti 



I. INTRODUCTION 

Subject to the effects of wind, a flat sandy surface is 
unstable and evolves into a regular periodic pattern of 
wavelength of the order of 10 cm and height of a few cen- 
timeters (see Fig. |l]) A slight sand accumulation 
on the surface tends to expose grains on the upwind stoss- 
side to the action of the wind, and shelter the grains on 
the downwind lee-side. Dislodged grains at the stoss tend 
to move toward the lee of the sand accumulation, the ini- 
tial perturbation is amplified, and small ripples develop, 
which migrate in the direction of the wind. Smaller rip- 
ples travel faster than larger ripples — since the migration 
velocity of ripples depends on the amount of grains be- 
ing transported during migration — so that small ripples 
merge with the larger ones. Due to ripple merging, the 
wavelength of the ripples grows in time as observed in 
the field [^| Jlo|l and in wind tunnel experiments |ll| , |l2|] . 

Due to the action of the wind, grains fly above the 
bed and strike the ground at small angles and with ve- 
locities given by the wind velocity. Successive impacts 
of grains are called saltation. As a result of bombard- 
ment of saltating grains, a number of ejected grains are 
generated from the static bed which subsequently move 
by surface creep. These grains (called restating grains) 
move distances smaller than the typical saltating jump 
(which is controlled by the wind velocity) until they are 
captured and become part of the bed. 

Climbing of ripples as seen in Fig. ^| is described by a 
vector of translation which has two components the 
component in the horizontal direction is the rate of ripple 
migration across the sediment surface, and the compo- 
nent in the vertical direction is the net rate of deposition, 
defined as the rate of displacement of the sand surface in 
the vertical direction (see Fig. ||a). The angle of climb- 
ing (£ in Fig. ||a) defines different types of lamination 
structures (see Chapter 9 0). If the angle of climb is 
smaller than the slope of the stoss, then the structures 
are called "subcritically" climbing ripples, and these are 



the sedimentary structures which are the focus of this 
study. 




FIG. 1. Wind sand ripples on a sandy surface in Death 
Valley National Park, California. 

Due to the climbing of ripples and segregation of grains 
during deposition and transport, ripple deposits develop 
lamination structures which are later preserved during 
solidification of the rock [||6[0]. The basic types of the 
smallest stratification structures in climbing ripples and 
small aeolian dunes relevant to this study are summa- 
rized by Hunter j| (see also J?J). In this paper we will 
focus on two lamination structures which are commonly 
found in sedimentary rocks: 

• Inverse-graded climbing ripple lamination: one of 
the most common lamination structures which are 
formed because grains composing the ripples dif- 
fer in size Jl|]^,|7|,|l3| . Large ejected grains travel in 
shorter trajectories than small grains, so that small 
grains are preferentially deposited in the trough of 
the ripples while large grains are deposited prefer- 
entially near the crest. Due to this segregation ef- 
fect, the migration of a single ripple produces two 
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layers of different grain size parallel to the climbing 
direction of the ripple as shown in Fig. ||a. This 
lamination structure is called inverse-graded — ac- 
cording to the nomenclature of Hunter ||] — since 
in a pocket of two consecutive layers formed by the 
climbing of a single ripple, the layer of large grains 
is above the layer of small grains (Fig. ||a). 

• Cross-stratification: Migration of ripples also pro- 
duces successive layers of fine and coarse grains 
not in the direction of climbing but parallel to 
the downwind face of the ripples (see Figs. || 
and 

1°) S!IH§- These structures— called cross- 
stratification or foresets — are also coarser toward 
the troughs ]l|,|i4|-[l6| as opposed to the segregation 
in inverse-graded climbing ripples. 




FIG. 2. Diagram showing the climbing of the ripples and 
the formation of cross-stratification patterns (from o). 

The instability giving rise to aeolian ripple morpholo- 
gies has been the subject of much study. The classic 
work of Bagnold jjj has been followed up by a number 
of studies. Modern models of wind ripple deposits are 
usually defined in terms of the splash functions intro- 
duced by Haff and coworkers [ [l7||l8|l (see also jl2| for a 
review) in order to model impact processes. These mod- 
els have successfully reproduced the instability leading 
to ripple deposits when the sand grains are of only one 
size 



When such models are defined for two 
species of grains differing in size [ p~3| | ? they result in pat- 
terns which resemble very closely the stratigraphic pat- 
terns found in inverse grading climbing ripple structures. 
It was in this context that Forrest and Haff proposed j^J 
that grading changes in ripple lamination are related to 
fluctuations in wind or transport. However, Anderson 
and Bunas Q found using a cellular automaton model 
that inverse-graded ripple lamination is due to different 
hopping lengths of small and large grains. These models 
do not include the interactions between the moving grains 
and the static surface j|(J — i.e., grains are assumed to 
stop as soon as they reach the sand surface — which are 
expected to be relevant for the dynamics in the rolling 
face of the ripples pl| , p2| . Recent studies by Terzidis et 
al |2y] have reproduced the ripple instability using the 



theory of surface flows of grains proposed by Bouchaud 
and collaborators This theory takes into account 

the interaction between rolling and static grains. More- 
over, recent studies of avalanche segregation of granular 
mixtures show dramatic effects when such interactions 
are taken into account |24U: 
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FIG. 3. (a) Cross-section of a sandstone showing 
small-scale lamination in 'subcritically' (as defined in [^| and 
[0]) inverse-graded ripples (along with the quantities defined 
in the model) . Each pair of layers of small and large grains is 
produced by the climbing of a single ripple, (b) Example of 
cross-stratification. Each climbing ripple produces a series of 
layers of small and large grains oriented across the direction 
of climbing, parallel to the downwind face. The basic types 
of the smallest stratification structures in ripples and small 
aeolian dunes relevant to this study are summarized in M. 



In this article, we use this formalism to include also 
the segregation effects arising when considering two type 
of species of different size and shape. We first formu- 
late a discrete model for two-dimensional transverse aeo- 
lian climbing ripples which incorporates simple rules for 
hopping and transport. Then, we incorporate the differ- 
ent properties of the grains, such as size and roughness, 
and we show that pattern formation in ripple deposits 
arises as a consequence of "grain segregation" during the 
flow and collisions. Specifically we show that three segre- 
gation mechanisms contribute to layer formation during 
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ripple migration (Fig. ||): 

• Size segregation due to different hopping lengths of 
small and large grains: the hopping length of the 
large reptating grains is smaller than the hopping 
length of small grains. 

• Size segregation during transport and rolling: 
larger grains tend to roll down to the bottom of 
troughs while small grains tend to be stopped pref- 
erentially near the crest of ripples. 

• Shape segregation during transport and rolling: 
rounded grains tend to roll down easier than more 
faceted or cubic grains, so that the more faceted 
grains tend to be at the crest of the ripples, and 
more spherical grains tend to be near the bottom. 

As a result of a competition between these segrega- 
tion mechanisms a richer variety of stratigraphic pat- 
terns emerges: inverse-graded ripple lamination occurs 
when segregation due to different jump lengths domi- 
nates, and cross-stratification (and normal-graded ripple 
lamination as well) occurs when size and shape segrega- 
tion during rolling dominates. Thus, a general framework 
is proposed to unify the mechanisms underlying the ori- 
gin of the most common wind ripple deposits. In the case 
of inverse-graded climbing ripple lamination we confirm 
a previously formulated hypothesis by Bagnold Q and 
Anderson and Bunas |l3| that the lamination is due to 
hopping induced segregation due to collisions. Regarding 
the origin of cross-stratification structures, we show that 
they arise under similar conditions as in recent table-top 
experiments of avalanche segregation of mixtures of large 
faceted grains and small rounded grains poured in a ver- 
tical Hele-Shaw cell ||,||. 

Segregation due to different hopping 
lengths 
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Size segregation: large grains at the bottom 
of the rolling-face and small grains at the top 



FIG. 4. Three segregation mechanism acting when the 
grains differ in size and shape. 



In what follows we first define the model for the case of 
single species of grains. We show that the model predicts 
the formation of a ripple structure, and propose a sim- 
plified continuum theory to understand the onset of the 
stability leading to ripple structure. Then, we generalize 
to the case of two type of species of grains differing in size 
and shape which give rise to the segregation mechanisms 
and the lamination structures seen in aeolian climbing 
ripples. 



II. ONE SPECIES MODEL 

We start by defining the model for the case of one 
type of grain in a two-dimensional lattice of lateral size 
L with periodic boundary conditions in the horizontal re- 
direction. Our main assumption is to consider two differ- 
ent phases [p2[-p6[ : a reptating or rolling phase composed 
by grains moving with velocity v by rolling, and a static 
phase composed by grains in the bulk. We also consider 
a curtain of external saltating grains which impact — at 
randomly chosen positions on the static surface — from 
the left to the right at a small angle f3 to the horizontal 
(Fig. |3|a). Shadow effects are considered by allowing only 
impacts with ballistic trajectories which do not intersect 
any prior portion of the surface. Upon impacting on the 
sand surface, a saltating grain dislodges n rep grains from 
the static surface, and jumps a distance l sa i after which 
the saltating grain is incorporated to the reptating phase. 
The n rep grains dislodged by the saltating grain jump a 
distance l rep <C l sa i- Upon reaching the surface, the dis- 
lodged grains form part of the reptating phase and move 
with velocity v. 

At every time step (At) only one reptating grain (in 
contact with the surface) interacts with the static grains 
of the bulk according to the angle of repose r p7| ; the 
remaining reptating grains move a distance vAt to the 
right. The angle of repose is the maximum angle at which 
a reptating grain is captured on the sand bed. If the lo- 
cal coarse-grained angle of the surface is smaller than the 
angle of repose, < r , the interacting reptating grain 
will stop and will be converted into a static grains. If the 
angle of the surface is larger than the angle of repose, 
> r , the reptating grain is not captured — and moves 
to the right a distance vAt with the remaining reptating 
grains — but ejects a static grain from the bulk into the 
reptating phase. 

The model predicts the formation of a ripple structure 
as seen in aeolian sand formations. The onset of the in- 
stability leading to ripples occurs with ripples of small 
wavelength. However, due to the fact that smaller rip- 
ples travel faster than larger ripples — since smaller rip- 
ples have less amount of material to transport — merging 
of small ripples on top of large ripples is observed. This 
gives rise to the change of the characteristic wavelength 
as a function of time, as is observed in wind tunnel exper- 
iments [p"T|]l2| , as well as in field observations [p|Jlo|] . We 
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find that our model predicts an initial power law growth 
of the wavelength of the ripples X(t) ~ t 0A (see Fig. ||). 
The wavelength of the ripples seems to saturate after this 
power law growth to a value determined by the saltating 
length. Moreover, we also observe a subsequent growth of 
the wavelength (not shown in the figure) up to a new sat- 
uration value close to a multiple of the saltation length. 
This process continues and a series of plateaus are ob- 
served until one large ripple (or dunes) of the size of the 
system is formed. The series of plateaus was also ob- 
served in a recent lattice model of ripple merging [ p8[ . 
However, we believe that the existence of plateaus might 
be due to an artifact of the discreteness of the lattice, 
and may have no physical meaning: the rules associated 
with the derivatives of the sand surface break down when 
the local slope is bigger than a given value determined 
by the discretization used to define the slope of the sur- 
face. A similar artifact was observed in a discrete model 
of ripple morphologies formed during surface erosion via 
ion-sputtering in amorphous materials |29|. Moreover, we 
will see that when we introduce two species of grains, the 
wavelength saturates at a constant value and the plateaus 
disappear. 



1(T 



OX) 

a 

> 



O Field measurements, Ref. [10] 

■ Wind tunnel experiments, Ref. [11] 

♦ Wind tunnel experiements, Ref. [12] 

— Present model 



10 
10 




10 10 10 

time (min) 



10 



FIG. 5. Growth of the wavelength of the ripples accord- 
ing to our model and comparison with field and experimental 
observations. The number of impacts in the model is trans- 
lated to time using a typical impact rate of 10 7 impacts m~ 2 
s _1 jl^]. The simulation data have been shifted vertically 
by a multiplicative factor, so that only the general trend (the 
slope) of the curve should be compared with the experimental 
data. 

In Fig. H we compare the prediction of our model with 
the available experimental data from field and wind tun- 
nel experiments. We observe a fair agreement between 
our model and the experiments. However, we believe that 
this agreement is not conclusive. In fact, the same exper- 
imental data can be also fitted with the same accuracy 
by a logarithmic growth of the wavelength of the ripples 
as shown recently by Werner and Gillespie (!(]] , who pro- 
posed a discrete stochastic model of ripple merging and 



found a logarithmic growth of the wavelength (see also 



A. Continuum formulation of the model 

The onset of the instability leading to the initial rip- 
ple structure can be studied using the continuum theory 
proposed by Bouchaud et al. [^3| to study avalanche in 
sandpiles. Recently, the set of coupled equations for sur- 
face flow of grains of has been adapted to the problem 
of the ripple instability Here, we propose a version 
of this theoretical formalism suitable to the physics of our 
model to study the initial ripple formation. Let R(x, t) 
describe the amount of reptating grains in the rolling 
phase and h(x, t) the height of the static bed at position 
x and time t, which is related to the angle of the sand 
surface by 6(x,t) = —dh/dx for small angles. Our set of 
convective-diffusion equations for the rolling and static 
phases are the following 123] : 



dR(x,t) 
dt 



dR 

dx 



D 



d 2 R 

dx 2 



T(R,6), 



dh(x, t) 
dt 



= -T(R,6), 



(la) 



(lb) 



where v is the drift velocity of the reptating grains along 
the positive x axis taken to be constant in space and time, 
and D a diffusion constant. The interaction term T takes 
into account the conversion of static grains into rolling 
grains, and vice versa. We propose the following form of 
r consistent with our model: 

T(R, 0) = a[f3 - 6{x, t)} + *y[0(x, t) - r }R(x, t), (2) 

where a and 7 are two constants with dimension of veloc- 
ity and frequency respectively, a is proportional to the 
number of collisions per unit time of the saltating grains 
with the sand bed, and also to the number of ejected 
grains per unit time, while 7 is the number of collisions 
per unit time between the reptating grains and the sand 
bed when they creep on the downwind side of the ripple. 
The first term in (g) takes into account the spontaneous 
(independent of R{x,t)) creation of reptating grains due 
to collisions of the saltating grains — a process which is 
most favorable at angles of the bed smaller than (3. The 
second term takes into account the interaction of the rep- 
tating grains with the bed of static grains: the rate of the 
interaction is proportional to the number of grains R(x, t) 
interacting with the sand surface. Rolling grains become 
part of the sand surface if the angle of the surface 6(x, t) 
is smaller than the repose angle 6 r ("capture"), while 
static grains become rolling grains if 0{x, t) is larger than 
r ("amplification"). Higher order terms are neglected 
since we are interested in the linear stability analysis de- 
termining the origin of the ripple instability. 
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Insight into the mechanism for ripple formation is ob- 
tained by studying the stability of the uniform solution 
ofEqs. (0): 



R = ap/Wr) 



h = 0. 



(3) 



We perform a stability analysis by looking for solu- 
tions of the type: R(x, t) = Ro + R(x, t) and 0(x, t) = 
60 + 6(x, t). The linearized equations for R and are 



dR{x,t) dR 
dt dx 



-7^-v m e- ~fO r R + D 



d 2 R 

dx 2 ' 



ggOM) _ _^ s# _ ^ as 



(4a) 



(4b) 



where the migration velocity of the traveling wave solu- 
tion is 



v m = a(l - p/0 r ), 



(5) 



which indicates that the ripple velocity is proportional to 
the number of collisions per unit unit of saltating grains 
with the sand bed, and to the number of ejected grains 
per collision. 

By Fourier analyzing Eqs. (||) we find a set of two 
homogeneous algebraic equations for R and with non- 
trivial solutions only when the determinant of the result- 
ing 2x2 matrix is zero. We obtain a dispersion relation 
w±(k) where the two branches ± correspond to the solu- 
tions of the resulting quadratic equation for w(k). Then 
we take the limit v m /v <C 1, which corresponds to the 
physical fact that the translation velocity (||) is much 
smaller than the rolling velocity of the individual grains 
(see [p2[), and we arrive to the following dispersion rela- 
tion for w_(k) (w+(k) gives rise only to stable modes): 



Im[w_(k)] 



-(VmM (jO r V 2 ) k 2 



( 7 6> r ) 2 + k 2 (2~/0 r D + v 2 ) + D 2 k 4 



■oC—) 2 

v 



The asymptotic forms for small and large k are 



Im[w_(k)] 
lm[w_(k)] 



VmV , 2 

D 2 k 2 ' 



k 
k 



(6) 



(7) 



oo 



which indicates that the branch w_(k) corresponds to 
only unstable modes w_(k) < 0. Similar stability anal- 
ysis was performed by Anderson [|o| using a continu- 
ity model neglecting the interaction between rolling and 
static grains and by Terzidis et al [^2[ with the contin- 
uum model Eqs. (|l|) but with a different interaction term 
than Terzidis et al. find a band of unstable modes 
until a given cut off wave vector at large k. This be- 
haviour is due to higher order derivatives of 9 appearing 
in the interaction term used in p^ |. The most unsta- 
ble mode k* in our model is given by dw_(k*)/dk = 0, 



with k* = y^fOr/D which gives an estimate of the initial 
wavelength of the ripples: 



(8) 



The final wavelength will be determined by higher or- 
der nonlinear terms arising from a more complicated in- 
teraction term than the one used in (||). As mentioned 
above, using the full discrete model we find that after the 
appearance of initial small undulations, the wavelength 
of the ripples grows due to ripple merging. In the discrete 
model, the wavelength seems to saturate to a finite value, 
although this value seems to be determined by the finite 
size of the simulation system (we use periodic boundary 
conditions in the horizontal direction). 



III. A MODEL FOR TWO SPECIES DIFFERING 
IN SIZE AND SHAPE 

Next we generalize the model to the case of two type of 
species differing in size and shape. The segregation mech- 
anisms discussed in the introduction are incorporated in 
the model as follows: 

• Size segregation due to different hopping lengths: 
we define l a (to replace l rep of the mono disperse 
case) as the distance a reptating grain of type a 
travels after being collided by a saltating grain. If 
we call the small grains type s and the large grains 
type Z, then li < l s . This effect was incorporated 
in the discrete stochastic model of Anderson and 
Bunas H using a generalization of the splash func- 
tions proposed by Haff Q . 

The interaction between the rolling grains and the sur- 
face is characterized by four different angles of repose; 
6 a/ 3 for the interaction of an a reptating grain and a (3 
static grain (replacing 6 r of the mono disperse case). The 
dynamics introduced by the angle of repose are relevant 
at the downwind face where two extra mechanisms for 
segregation of grains act in the system: 

• Size segregation during transport and rolling: large 
rolling grains are found to rise to the top of the rep- 
tating phase while the small grains sink downward 
through the gaps left by the motion of larger grains; 
an effect known as percolation or kinematic sieving 
[0j3O|. Due to this effect only the small grains in- 
teract with the surface when large grains are also 
present in the rolling phase. Small grains are cap- 
tured first near the crest of the ripples causing the 
larger grains to be convected further to the bottom 
of the ripples. We incorporate this dynamical seg- 
regation effect by considering that, when large and 
small grains are present in the reptating phase, only 
the small ones interact with the surface according 
to the angle of repose, while the large grains, be- 
ing at the top of the reptating phase do not interact 



5 



with the static grains and they are convected down- 
ward. Thus, the large grains interact with the sur- 
face only when there are no small grains present in 
the reptating phase. A similar percolation mecha- 
nism was introduced in the discrete and continuum 
models of (24| to understand the origin of strat- 
ification patterns in two-dimensional sandpiles of 
granular mixtures. 

• Shape segregation during transport and rolling: 
rounded grains tend to roll down easier than more 
faceted or cubic grains, so that the more faceted 
grains tend to be at the crest of the ripples. This 
segregation effect is quantified by the angles of re- 
pose of the pure species, since the repose angle is 
determined by the shape and surface properties of 
the grains and not by their size: the rougher or the 
more faceted the surface of the grains the larger the 
angle of repose. If the large grains are more faceted 
than the small grains we have SS < On , while when 
the small grains are more faceted On < SS J24|] . The 
species with the larger angle of repose have a larger 
probability to be captured at the sand bed than the 
species with smaller angle of repose. 

We notice that the fact that the grains have different 
size leads to different cross-angles of repose 6 a p |^ , ^T| . 
If the large grain are type I and the small are type s then 
6is < @sh which in turn leads to another size segregation 
effect. However, by incorporating the size segregation 
due to the different cross-angle of repose we would get 
the same effect as the one incorporated in the model by 
the percolation effect (the small grains are preferentially 
trapped at the top of the slip-face), so that, in what 
follows, we will consider only the percolation effect for 
simplicity — i.e., large reptating grains are allowed to in- 
teract with the sand surface only when there are no small 
reptating grains below the large grains at a given posi- 
tion. 

We also notice that, in general, the distance a grain 
travels after being kicked by a saltating grain depends 
on the type of colliding grain and the type of grain on 
the bed. Thus we define l a p as the distance a reptating 
grain of type (3 travel after being collided by a saltating 
grain of type <x If we call the large grains type I and the 
small grains type s then we have: l s i < In < l ss < li s . 
Hoverer, in practice we will make the following approx- 
imation I s i = In = h, an d lss — ks = i-e., we do not 
consider which type of grain is colliding. 

Thus, three different mechanisms — size-segregation 
due to percolation in the reptating phase, shape- 
segregation in avalanches, and size segregation due to 
the different reptating jumps l a — compete in the system 
giving rise to a rich variety of lamination patterns as we 
show below. 



A. Size segregation due to different hopping lengths 

We first investigate the morphologies predicted by the 
model in the case where the different reptating jumps l a 
play the dominant role in the segregation process. This 
is the case l s — h > h and when the grains have ap- 
proximately the same shape SS « On. We start our sim- 
ulations from a flat sand surface composed by a 50/50 
by volume mixture of small and large grains and we ob- 
serve the system to evolve into ripples traveling in the 
direction of the wind (Fig. ||). The resulting morphol- 
ogy seen in Fig. || resembles the most common climbing 
ripple structures such as those documented in field ob- 
servations [Hj6[j|] (Fig. ||a). This result confirms the 
hypothesis of Anderson and Bunas j^J that the origin 
of inverse-graded lamination in climbing ripples is due to 
the size segregation effect produced by the different hop- 
ping lengths of small and large grains. In our simulations 
we observe that large grains (dark color) are deposited 
at the top of the ripples while small grains (light color) 
are deposited preferentially at the bottom of the ripples 
since k < l s , resulting in a lamination structure showing 
inverse grading (layer of large grains on top of the layer 
of small grains). 

The system size used in our simulations corresponds 
to 256 bins, with At = 1 s, n rep = 4 grains per impact, 
I sal = 90 cm, (3 = 10°. For the case shown in Fig. |6| we 
use l\ = 25 cm and I2 = 5 cm, and Q a p = 30°. 

Inverse graded lamination 



wind 




FIG. 6. Morphology predicted by the model when segre- 
gation due to different jump lengths is dominant showing in- 
verse grading climbing ripple lamination. From left to right, 
we show a sequence of three stages in the dynamics of the 
climbing ripples. Starting from a flat sand surface of small 
(light color) and large grains (dark color), the system evolves 
into ripples climbing in the direction of the wind. 

The efficiency of the segregation mechanisms is sensi- 
tive to the relative value of the parameters of the model. 
For instance we find that by decreasing the value of l sa i 
relative to the reptating lengths, the efficiency of segrega- 
tion is greatly reduced in the case of inverse-graded lam- 
ination. The structures shown here are all subcritically 
climbing from left to right (for a review of the termi- 
nology and definition of different climbing structures in 
ripples deposit see ||], and Chapter 9 0]). Supercritical 
climbing ripples are much less common and are produced 
in slowly translating ripples with small number of ejected 
grains per impact |l9| ]. We also notice that the shape of 
ripples found in Nature are more asymmetric with the 
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downwind side at a steeper angle than the upwind side, 
while our model predicts a more symmetric triangular 
shape of the ripples. More realistic asymmetric shape of 
the ripples can be obtained by considering the exact tra- 
jectories of the flying grains and the complex interaction 
of grains with the air flow as shown by previous models 
by Anderson Q. 



B. Segregation due to rolling and transport 

Next we consider the case where the difference between 
the reptating jumps is small, \l s — k\/l s < 1 (or when 
the downwind face is large compared to l s or l{) and the 
grains differ appreciable in shape (0 SS ^ On) and also 
in size. Then segregation in the rolling face is the rele- 
vant mechanism for segregation and we do not take into 
account the segregation due to different hoping lengths. 

We first consider the case where the small grains are 
the roughest and the large grains are the most rounded 
(0 SS > On). I n this case a segregation solution along the 
downwind face is possible since both, size segregation due 
to percolation and shape segregation act to segregate the 
small-faceted grains at the top of the crest and the large- 
rounded grains at the bottom of troughs. The result is a 
lamination structure (Fig. |?]a) which resembles the struc- 
ture of climbing ripples of Fig. || but with the opposite 
grading: small grains at the top of each pair of layers 
(normal-grading climbing ripple lamination). This type 
of lamination is not very common in Nature. 

On the other hand, when the large grains are rougher 
than the small grains (On > SS ), a competition be- 
tween size and shape segregation occurs Q. Size seg- 
regation due to percolation tends to segregate the large- 
cubic grains at the bottom of the downwind face, while 
the shape segregation effect tends to segregate the same 
grains at the crest of the ripples. Then a segregation 
solution along the downwind face as in Fig. 0a is not 
possible, and the result of this instability is the appear- 
ance of layers of small and large grains parallel to the 
downwind face (Fig. ^|b) and not perpendicular as in 
Fig. ^a. These structures correspond to the more com- 
mon cross-stratification patterns found in rocks 
In addition to the stratification parallel to the rolling 
face, the deposits are also coarser toward the bottom of 
the downwind face. 

The mechanism for cross-stratification involves two 
phenomena superimposed, (i) Segregation in the rolling 
face: a pair of layers is laid down in a single rolling event, 
with the small grains segregating themselves underneath 
the larger grains through a "kink" mechanism as seen in 
]l6| , p4[ . This kink in the local profile of static grains pro- 
vides local stability to the rolling grains trapping them. 
The kink moves uphill forming layers of small and large 
grains, (ii) The climbing of ripples due to grain deposi- 
tion and ripple migration. 

The conditions for cross-stratification and the dynam- 



ical segregation process found with the present model 
are similar to the findings of the experiments of strati- 
fication of granular mixtures in two-dimensional vertical 
cells performed in @ and also the models of (24) . We 
observe the formation on the slip face of the ripples of 
an upward traveling wave or "kink" at which grains are 
stopped during the avalanche as was observed in ]l6|,Q . 
According to this, the wavelength of the layers is pro- 
portional to the flux of grains reaching the rolling face, 
then it is proportional to the number of saltating grains 
impacting the surface. 

For the case shown in Fig. [?]a we use l 3 = = 5 cm, 
and s i = SS = 30°, and is = M = 20°, while for Fig. 
0b we use l s = k = 5 cm, and s i = SS = 20°, and 
ls = U = 30°. 

(a) Normal graded lamination 




(b) Cross-stratification 




FIG. 7. Resulting morphologies predicted by the model 
after 10 7 impacts, when segregation in the rolling face is 
dominant showing (a) normal grading lamination of large 
rounded grains (dark color, at the bottom of the downwind 
face) and small rough grains (light color, at the top); and 
(b) cross- stratification of large rough grains (dark color) and 
small rounded grains (light color). 
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C. Phase diagram: General case 



We have also investigated the morphologies obtained 
when the three segregation effects (shape segregation, 
percolation effect, and hopping induced size segrega- 
tion) act simultaneously. The resulting morphologies are 
shown in Fig. along with the phase diagram summa- 
rizing the results obtained with the model in Figs. and 

0- 

The case \l 3 — k\ > l 3 and SS « On corresponds to the 
inverse-graded lamination shown in Fig. 0, as indicated 
in the phase diagram Fig. |[ Moreover, when l s — h > l s 
and for any other value of the angles SS and On , we find 
that the hopping induced segregation seems to dominate 
over rolling induced segregation. As can be seen from the 
upper two panels in Fig. 0, the morphologies obtained 
in these cases resemble the one of Fig. corresponding 
to inverse-graded climbing ripple lamination. The upper 
left panel shows the results when the smaller grains are 
the roughest On — S s< and shows clearly the same 
lamination as in Fig. |6[ In this case the hopping induced 
segregation dominates completely over the rolling segre- 
gation due to the angles of repose and percolation. In 
the upper right panel we show the case when the large 
grains are the roughest On — SS > 0. In this case we 
also see the same lamination pattern as in Fig. © but we 
see a tenuous trace of the cross-stratification seen in Fig. 
0b too; the rolling induced segregation seems to have a 
more important role in comparison with the hopping in- 
duced segregation than in the case shown in the upper 
left panel. 

The region \l s — li\ < l s is discussed in Figs. 0, 
and our model shows normal grading lamination when 
On — SS < (Fig. 0a) and cross-stratification when 
M -0 M >O(Fig. 0b). 

Interesting morphologies are predicted when l 3 — h < 
—l s . This might be the case of a large difference in density 
between the grains, i.e., very light large grains and heavy 
small grains. The patterns obtained with the model are 
shown in the two lower panels in Fig. 0. The late stage 
in the dynamics evolution of the patterns shown at the 
left lower panel (On — SS < 0) and at the right lower 
panel (On — SS > 0) in Fig. resemble the patterns 
of normal grading climbing ripples (Fig. 0a) and cross- 
stratification (Fig. 0b) respectively. This indicates the 
dominance of avalanche segregation (shape segregation 
and percolation effect) over the hopping induced segre- 
gation in this region of the phase space (this dominance 
is the opposite of what we find in the region l s — li > l s ). 
However, as seen in Fig. 0, these patterns appear only 
in the late stages of the evolution. As seen in the two 
lower panels, in the early stages the lamination structure 
appears to be dominated by the hopping induced seg- 
regation. We observe a clear transition at intermediate 
time in the simulation which can be clearly observed in 
the change of the climbing angle of the ripples observed 
in the left and right lower panels in Fig. p| The late stage 



patterns appear only when the slip-face is well developed 
after a transient of small ripples dominated by hopping 
segregation. We are not aware of any wind tunnel ex- 
periment or field observation showing similar dynamical 
transition, so it would be interesting to explore this tran- 
sition experimentally. 




normal grading 



climbing ripple, Fig. 7a 




cross-stratification 



Fig. 7b 




FIG. 8. Phase diagram predicted by the model. The "s" 
refers to the small grains (light colors), and the "/" refers to 
the large grains (dark colors). SS > On means that the smaller 
grains are the roughest, and On > SS means that the larger 
grains the roughest. The three areas in the phase space near 
the axis refer to the morphologies already studied in Figs. ^, 
0a, and ^b. 



IV. SUMMARY 

In summary, we have shown that lamination in sedi- 
mentary rocks at around the centimeter scale is a mani- 
festation of grain size and grain shape segregation. Thin 
layers of coarse and fine sand are present in these rocks, 
and understanding how layers in sandstone are created 
might aid, for instance, in oil exploration since great 
amounts of oil are locked beneath layered rocks. Our 
findings suggest a unifying framework towards the under- 
standing of the origin of inverse grading, normal grading 
and cross-stratification patterns in wind ripple deposits. 
We identify the conditions under which these different 
laminations patterns arise in sandstone. In particular, 
we find that cross-stratification is only possible when the 
large grains are coarser than the small grains composing 
the layered structure in agreement with recent experi- 
ments on avalanche segregation in two-dimensional sand- 
piles. The fact that cross-stratification patterns are very 
common in Nature might be due to the fact that frag- 



8 



mentation processes usually leads to smaller grains more 
rounded than larger grains. So far we have explored a 
two-dimensional cross-section of the ripple deposits. Ex- 
tensions to three dimensional systems are ongoing as new 
physics may emerge when taking into account the lateral 
motion of the grains |2l| . 
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